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ABSTRACT 

We use 3D radiative MHD simulations to investigate the formation and dy- 
r£3 ■ namics of small-scale (less than 0.5 Mm in diameter) vortex tubes spontaneously 

generated by turbulent convection in quiet-Sun regions with initially weak mean 
magnetic fields. The results show that the vortex tubes penetrate into the chro- 
mosphere and substantially affect the structure and dynamics of the solar at- 
mosphere. The vortex tubes are mostly concentrated in intergranular lanes and 
are characterized by strong (near sonic) downflows and swirling motions that 
capture and twist magnetic field lines, forming magnetic flux tubes that expand 
with height and which attain magnetic field strengths ranging from 200 G in the 
chromosphere to more than 1 kG in the photosphere. We investigate in detail the 
physical properties of these vortex tubes, including thermodynamic properties, 
flow dynamics, and kinetic and current helicities, and conclude that magnetized 
vortex tubes provide an important path for energy and momentum transfer from 
the convection zone into the chromosphere. 

Subject headings: Sun: photosphere, chromosphere, surface magnetism, magnetic 
topology 



1. Introduction 



Interest in vortex tube dynamics of the quiet Sun was recently i nitiated by the de 
tection of ubiquitous small-scale swirling motions in the photosphere (jWang et al. 1 11995 
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Potzi fc Brandt 2005 



Bonet et al. 


2008, 


2010; 


Balmaceda et al. 


2010; 


Steiner et al. 


2010) 


(Wedemever-Bohm & Rouppe van der Voort 


200< 


)) with high-resoluti 



solar telescopes. Pr evious to this discovery, vortex tubes on the S un were predicted by theo- 



retical models (e.g. , IStenflo 1119751) and numerical simulations (e.g-. lBrandenburg et al. 111996 



Stein fc Nordlund 200oh . giving a clear illustration of the turbulent nature of solar convec- 
tion. Both observations and numerical simulations show concentrations of vortex tubes in 
the intergranular lanes. According to rec ent radiative hydro-dyna mic simulations, vortical 
motions can be also form inside granules ( jKitiashvili et al. 1 120121 ) . These simulations have 
also shown that vortex tube formation in the near-surface layers can be caused by two basic 
mechanisms associated with: 1) small-scale convective instability developing inside granules, 
and 2) the Kelvin-Helmholtz instability of shearing flows. 

The convective instability leads to formation of a vortex sheet and its subsequent over- 
turning during a localized upflow (plume) or splitting of a granule. The process of the 
vortex sheet overturning, which results in a vortex tube, is often accom panied by a grad- 



ual m igration of the vortex tube into an intergranular lane (see Fig. 2 in iKitiashvili et al. 



20121 ). Shearing flows that lead to the development of the Kelvin-Helmholtz instability can 
be present in both granules and intergranular lanes. However, in the intergranular lanes the 
shearing flows are stronger and can lead to a series of vortices (resembling the Karman vortex 
street). Also, converging downflows in the intergranular lanes make the vortex tubes more 
stable, with characteristic lifetimes up to 40 min, whereas inside granules the lifetime is less 
than 10 min. These processes can explain why the observed vortex tubes are predominantly 
concentrated in the intergranular lanes. 

Numerical simulations also show connections between vort ex tube dynamics and var i- 
ous other solar phenomena, such as the hydromagnetic dynamo ([Brandenburg et al. Ill996l ). 
spontaneous organizati on of emerged magnetic field into self-maintai ned pore-like structure s 
( jKitiashvili et al. excitation of acoustic waves in the quiet Sun (jKitiashvili et al.ll201ll ). 

and others. In this Letter, we present new numerical simulations that demonstrate important 
links between the turbulent subsurface layers and the solar atmosphere though the dynamics 
of penetrating vortex tubes. 



2. Computational setup 



Numerical simulations of the quiet Sun are performed by using a 3D radiative MHD 
code ('SolarBox') developed at the NASA/Ames Research Cent er and the Stanford Cen- 
ter for Turbulence Research by Alan Wray and his colleagues (IJacoutot et al.l 120081 ) for 
modeling the outer part of the solar convection zone and lower atmosphere in a carte- 
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sian geometr y. The code was de v eloped for realistic-type numerical simulations of the Sun 
pioneered by iNordlund fe Stein I (120011 ) and uses a tabular real-gas equation of stat. Ra- 



diative energy transfer is calculated with a 3D multi-spectral-bin method between fluid 



(Rogers et al. 


1996 


) . Initialization 


( Christensen-Dalsgaard et al. 


1996) 



The physical description of the dynamical properties of solar convection was improved 
through the implementation of subgrid-scale turbulence models, which effectively increase 
the Reynolds number and allow better resolution of essential turbulent scales. This ap- 
proach, based on Large-Eddy Simulation (LES) models of subgrid turbulence, has demon- 
strated good agreeme nt of numerically modeled acoustic wave excitation with observations 



( Jacoutot et al. 20081) and ha s helped improve understanding of wave excitation mecha - 



nisms (jKitiashvili et al.l 1201 ll ) , formation of magnetic structures ( jKitiashvili et al.l |2010| ) 



and Ever shed flows in sunspots (jKitiashvili et al.l 120091). The simu 



at ions in this paper were 



1963) in which the com- 



1991; Jacoutot et al. 



obtained using a Smagorinsky eddy-viscosity model (ISmagorinskv 
press ible Reynolds stresses were calculated in the form (IMoin et al. 

2008h : Tij = -2CsA 2 \S\(Sij - w fc , fc <%/3) + 2C c A 2 \S\ 2 S ij /3, where the Smagorinsky coeffi- 
cients Cs = Cc = 0.001, Sij is the large-scale stress tensor, and A = (dx x dy x dz) 1 ^ with 
dx, dy and dz being the grid-cell dimensions. 

In the current study, the simulation results were obtained for a computational domain 
of 6.4 x 6.4 x 6.2 Mm 3 , including a 1 Mm high layer of the atmosphere, with a grid spacing 
of dx = dy = 12.5 km and dz = 10 km. The lateral boundary conditions are periodic. The 
top boundary is open to mass, momentum, and energy transfers and also to radiative flux. 
The bottom boundary is open for radiation and flows, and simulates energy input from the 
interior of the Sun. We focus mostly on a case with an initially uniform vertical magnetic 
field, B z = 10 G, representing quiet-Sun conditions (far from sunspots and active regions). 



3. Formation of vortex tubes by turbulent convection 

Vorte x tubes are formed by turbulent convection in near -surface layers of the convective 



zone (e.g. IStein fc Nordlund 1 12000| ; iKitiashvili et al. 1 120121 ) . The vortex tubes represent 



compact low-density structures up to 0.5 Mm in diameter and with high-speed swirling 
motion reaching up to 12 km/s. The vortex cores are characterized by strong downflows 
(up to 8 km/s) and lower temperature. Large vortex tubes can extend deeper than 300 km 
below the surface. 
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Our previous simulations (iKitiashvili et al. 1 12012| ) revealed two basic mechanisms of 
vortex tube formation: one due to a granular instability (vortex sheet overturning) and 
another due to the Kelvin-Helmholtz instability in shearing flows. Vortex tubes can form in 
intergranular lanes and in granules but are mostly concentrated in the intergranular lanes 
(Fig. [TJa, b). These physical mechanisms of vortex tube formation are purely hydrodynamic, 
but in the real Sun vortices are expected to strongly interact with ubiquitous magnetic fields. 
However, neither simulations nor observations have shown a clear correlation between vortex 
motions and magnetic field concentrations, that is, not every vortex is accompanied by a 
strong magnetic field concentration. Th is fact has also been shown in simul ations using a 
shallow domain (1.4 Mm in total height; Ishelvag et al.ll201ll ; IMoII et al.ll201l[ ). 



In weak magnetic field regions, magnetic patches follow convective motions. Concen- 
tration and magnification of magnetic field by swirling motions can stabilize the vortex tube 
structure and decrease the influence of surrounding turbulent flows. In our simulation case, 
we introduce a 10 G, initially uniform, vertical magnetic field. This field gets quickly concen- 
trated, mostly in intergranular lanes, and we find that the strongest magnetic field (~ 1 kG) 
concentrations are often associated with vortices (Fig. [T]c, d). 



4. Dynamics and properties of vortex tubes in the chromosphere 

A new interesting result of our simulations is the extension of turbulent vortex tubes 
from the convection zone into the convectively stable atmospheric layers. Figure [2] illustrates 
a snapshot of enstrophy distribution showing vortex tube structures (yellow isosurfaces) 
above the photosphere (the horizontal wavy light surface shows the 6400 K near-surface 
layer). These vortex tubes are mostly concentrated in the intergranular lanes and often form 
arc-shaped structures above the surface. Other vortices penetrate almost vertically into the 
higher chromospheric layers (an example of such an extended vortex tube is indicated by the 
arrow; we will consider its structure in detail below). Local upflows (red color on vertical 
slices, Figure [2]) cause stretching of the vortex arcs, and nearby vortices can destroy them. 
Finally, propagating shock waves interact with the vortex tubes in the higher chromospheric 
layers. The overall chromospheric dynamics driven by turbulent convection is thus very 
complicated. The effect of vortex penetration into the chromosphere is mostly hydrodynamic, 
as observed in simulations with and without magnetic field. However, the magnetic field 
tends to be captured and concentrated in the vortex tubes, causing new dynamical effects. 

The structure of the vortex tube indicated by the arrow in Figure [2] is illustrated in 
Figure [3] at different heights: 200 km, 500 km, 650 km and 800 km above the surface. The 
temperature distribution (Fig. [31 row a) shows local heating of the vortex core region, whereas 
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in the subphotospheric layers the core vortex temperature is lower than in the surrounding 
plasma. For this moment of time, the swirling motions in the vortex region are characterized 
mostly by highly turbulent downflows (Fig. but some upflows are noticeable near the 
edge of the vortex tube, the size of which is expanding with height. The current helicity, 
calculated in Alfven units as Xm = ■ (V X B), forms a sheet-like structure oriented along 
the intergranular lane near the photospheric layers (Fig. [3](f). The current sheet structure 
gradually changes orientation in the higher layers and becomes more circular (Fig. |3]o?). The 
current helicity structure is more diffuse than the kinetic helicity, \ m — u- (V x u), shown in 
Figure [3]c. The density distribution in the lower atmosphere is similar to the surface layers, 
but the vortex tube structure becomes more complicated with height, forming a ring-like 
structure at ~ 800 km above the surface (Fig. [3]e). 

Figure H] shows the time evolution of the velocity streamlines (panels a-c), magnetic 
field lines (panels d-f), and the ratio of gas pressure to magnetic pressure (plasma /3) for 
three moments separated by 3 min. The structure of the vortex tube in the middle column 
is shown in more detail in Fig. [5]a. In Figure HJ the grey-yellow isosurface corresponds to a 
temperature of 5800 K. Color patches on this surface indicate variations of magnetic field 
strength as indicated in the right color bar. The strongest magnetic field concentrations 
(~ 1.2 kG) are associated with the vortex tubes in the photospheric layer; and the field 
strength decreases to ~ 200 G in the upper layers of our domain (~ 1 Mm above the 
photosphere). 

The numerical simulations show the penetration and dynamics of the vortex tube into 
the chromosphere. The vortex core contains very compact helical downflows, and we observe 
that the vortex pulls granular fluid upward which then reverses into the downflows (Fig. |4]a). 
The magnetic field at this stage of vortex evolution continues to concentrate in the vicinity 
of the vortex by following the swirling turbulent motions (Fig. H]c?). These strong helical 
flows capture and twist the magnetic field lines. Also, the helical magnetic loops formed 
by vortex tubes have a tendency to move upward due to local upflows near the vortex 
core. Three minutes later, the helical downflows have become more compact and stronger 
(Fig. H]6), but the vortex is affecting a larger surrounding area. We begin to see evidence of 
vortex decay when this vortex starts interacting with others by sharing with them a part 
of the downflow (Fig. Finally, during the next three minutes, the photospheric and 

chromospheric parts of the vortex tube become disconnected but still continue to evolve. 
Figure H]c shows remnants of the initially strong helical flows in the atmosphere. At this 
moment, they are still weakly helical and become captured by another growing vortex. The 
magnetic field lines also keep their helical topology and start to diffuse (Fig. @J. We show 
the plasma parameter /3 = ^ = 3 level as blue isosurfaces in Figures IHg-z; the value /3 — 1 is 
reached only in a small region of the vortex core. This parameter shows that magnetic effects 
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play a significant role in the photospheric layers of the vortex tube and that the region of 
their influence rapidly expands with height. At the decay stage of the vortex tube (Fig. HJi), 
magnetic effects are significant only in the upper layers. 

The relative role of kinematic and magnetic effects of the swirling motions is illustrated 
by the kinetic, Xk, and magnetic, Xm, helicities. An example of relative distribution for 
both helicities is shown in Figure |5]a (the kinetic helicity is in blue, and the current one 
is in pink) for the vortex tube that is indicated by the arrow in Figure HJ Blue and pink 
isosurfaces correspond to helicity values of —5000 cm/s 2 ; the current helicity is calculated in 
Alfven units and has the same dimension as the kinetic helicity. In Figure [5^, we also plot 
the temperature isosurface for 5800 K, which has a very compact structure of a complicated 
chiralical shape, expanding into the higher layers of the atmosphere. The distribution of the 
kinetic helicity is more compact than the current helicity, meaning that the swirling flows in 
the vortex tube are more compact than the twisted magnetic field lines. 

In general, the dynamics of the subsurface and near-surface layers is dominated by tur- 
bulent convective motions, while magnetic effects are noticeable in the small-scale magnetic 
flux concentrations in the intergranular lanes (magnetic flux tubes). In higher atmospheric 
layers, magnetic effects are stronger because of the fast decrease of gas pressure, which leads 
to expansion of the magnetic flux tubes. 

To investigate the properties of the magnetized vortex tube with height, we selected 
a region inside the —5000 cm/s 2 isosurface of the current helicity. In this region we have 
plotted the mean values of temperature, vertical and horizontal velocities, magnetic field, 
and kinetic and current helicities as a function of height for different moments of time 
with a cadence of 20 sec (Fig. &)-f). The values of temperature and density are shown as 
perturbations from and normalized by the mean values: (T vortex — T)/T and (p vor tex — p) /p. 
The temperature distribution shows a deficit in the convective layers of the vortex tube. 
Above the photosphere, the temperature in the vortex tube increases, and we can see heating 
in the vortex core (Fig. \5$>). Occasional temperature decreases above ~ 500 km reflect the 
dynamically oscillatory behavior of the vortex tube. The density distribution (Fig. |5]) shows 
an increase below the surface due to mass concentration around the vortex core, which has 
significantly lower density. Above the surface, the mean density perturbation in the tube 
first decreases and then increases above 200 km (Fig. [5]c). 

The vertical distribution of the mean velocity inside the vortex tube shows very different 
properties for the horizontal and vertical components. The mean horizontal speed (blue 
curves, Fig[5]tf) is almost constant along the vortex tube, with relatively small fluctuations in 
time around the mean speed of ~ 3.5 km/s. It is interesting that the mean horizontal speed, 
averaged over time and over the whole domain, (thick dark blue curve) shows a decrease 
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at 400 — 800 km above the surface, but inside the vortex tube there is no such decrease. 
In contrast to the horizontal speed, the vertical velocity component (red curves, Fig, |5]tf) is 
very dynamic and is characterized by predominant downflows; however, local upflows can be 
detected inside the vortex tube. 

The vertical component of magnetic field is significantly stronger than the magnitude 
of the horizontal field (Fig. [5]e). Both vertical (red curves) and horizontal (blue curves) 
fields show a similar tendency to decrease in the atmospheric layers, which is reflected in 
the expanding topology of the flux tube. At a height of about 500 km, the magnitude of 
the vertical component of magnetic field is smaller than the horizontal component because 
the magnetic field lines become more twisted by the vortex. The mean kinetic helicity (blue 
curves, Fig. \5jf) is significantly greater than the mean current helicity (red curves) because 
swirling motions in the tube are accompanied by strong downflows, while the strongest 
magnetic field is only weakly twisted. 

5. Conclusion 

The formation and dynamics of small-scale vortex tubes play key roles in various pro- 
cesses in solar surface convection and in the solar atmosphere. Our radiative MHD simu- 
lations reveal vortex tubes formed by turbulent convection penetrating from the subphoto- 
sphere into the chromosphere. These vortex tubes cause significant qualitative changes in 
atmospheric dynamics, leading to strong variations in the thermodynamic structure through 
local heating and density variations, generating twisted magnetic flux tubes, and creating 
local twisted upflows into the chromosphere. Strong localized swirling motions occupy large 
areas around the vortex tubes, capturing and twisting magnetic field lines from nearby 
magnetic structures. As a result of these phenomena, magnetized vortex tubes generated 
by turbulent convective motions provide a very important link for energy and momentum 
exchange between the surface layers and the chromosphere. 

This work was partially supported by the NASA grant NNX10AC55G, the International 
Space Science Institute (Bern) and Nordita (Stockholm). The authors thank Phil Goode, 
Vasyl Yurchshin, Valentyna Abramenko, and participants of the Nordita and ISSI teams for 
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Fig. 1. — Snapshots of vertical velocity V z and enstrophy u 2 = (V x u) 2 on the solar surface 
for two simulation cases: without magnetic field (panels a and b) and with initial uniform 
weak vertical magnetic field, Bz = 10 G (panels c and d). Examples of small-scale vortex 
tubes are indicated by white and black squares. 
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Fig. 2. — 3D snapshot of the top part of the computational domain shows the penetration of 
convective vortex tubes (yellow isosurfaces) from the subphotosphere into the low and mid 
chromospheric layers. The horizontal wavy surface indicates the distribution of temperature 
at 6400 K, and corresponds to a photosphere layer. The vortex tubes (yellow isosurfaces) 
are shown for the enstrophy value of 0.0075 s~ 2 . The vertical slices in the back illustrate 
the vertical velocity distribution. Blue color corresponds to downflows, yellow-red shows 
upflows. 
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Fig. 3. — Vortex tube structure in four different atmospheric layers: 200 km, 500 km, 650 km 
and 800 km above the solar surface (from left to right). Each row corresponds to various 
quantities: a) temperature T, b) vertical velocity Vz, c) kinetic helicity Xk, d) current helicity 
Xm, and e) density p. 



13 




Fig. 4. — Evolution of the velocity field (streamlines in panels a-c), magnetic topology 
(streamlines in panels d-e), and plasma parameter (3 (blue isosurface for (3 = 3 in panels 
g-i). Each column corresponds to simulation data 3 min apart. The grey isosurface shows 
T=5800 K; additional coloring from light yellow to orange indicates variations of the mag- 
netic field strength in the range from to 1200 G. Coloring of the velocity streamlines in 
panels a-c corresponds to vertical velocities from —7 km/s (blue) to +7 km/s (red). 
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Fig. 5. — a) 3D rendering of a vortex tube indicated by the arrow in Figure [2] shows the 
relative distribution of the kinetic (blue isosurface) and current helicities (pink) for value 
—5000 cm/s 2 . The grey- yellow isosurface shows the distribution of temperature T=5800 K. 
Panels b-f show the distribution with height of mean vortex tube properties at different 
moments of time with 20 sec cadence: b) relative temperature variations, c) relative density, 
d) horizontal speed (Vh) and vertical velocity (Vz), e) horizontal (Bh) and vertical (Bz) 
magnetic field strength, f) kinetic (xk) and magnetic (xm) helicities. Thick blue solid curve 
in panel d shows the mean horizontal speed averaged over the whole domain and time. The 
systematic variations at z > 0.95 Mm are due to the top boundary conditions. 



